Preprint typeset in JHEP style - HYPER VERSION 



On measuring the masses of pair-produced 
semi-invisibly decaying particles at hadron colliders 



Daniel R. Tovey 

Department of Physics and Astronomy, 
University of Sheffield, 
Hounsfield Road, Sheffield S3 IRE, UK 
E-mail. 



daniel . toveyOcern . ch 



Abstract: A straightforward new technique is introduced which enables measurement at 
hadron colliders of an analytical combination of the masses of pair-produced semi-invisibly 
decaying particles and their invisible decay products. The new technique makes use of the 
invariance under contra-linear Lorentz boosts of a simple combination of the transverse 
momentum components of the aggregate visible products of each decay chain. In the 
general case where the invariant masses of the visible decay products are non-zero it is 
shown that in principle the masses of both the initial particles from the hard scattering and 
the invisible particles produced in the decay chains can be determined independently. This 
application is likely to be difficult to realise in practice however due to the contamination 
of the final state with ISR jets. The technique may be of most use for measurements of 
SUSY particle masses at the LHC, however the technique should be applicable to any class 
of hadron collider events in which heavy particles of unknown mass are pair-produced and 
decay to semi-invisible final states. 
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1. Introduction 

In R-Parity conserving SUSY events at hadron colliders SUSY particles ('sparticles') must 
be pair-produced and undergo cascade decay to the Lightest Supersymmetric Particle 
(LSP), which is often invisible and hence a dark matter candidate. The presence of two 
such invisible particles in the final state, together with imperfect detector hermeticity close 
to the beam-pipe and an uncertain parton centre-of-mass energy, prevents the use of con- 
ventional invariant mass or transverse mass techniques for sparticle mass measurement. 
Similar challenges are faced when attempting to measure the mass of any pair-produced 
particles with visible and invisible decay products. 

Several approaches to this general problem have been documented, usually in the con- 
text of measuring SUSY particle masses. Given a sufficiently long decay chain constraints 
on analytical combinations of sparticle masses can be obtained from the positions of end- 
points in distributions of invariant masses of combinations of visible SUSY decay products 
(jets, leptons etc.) Q. Given a number of such constraints, the system of equations may be 
solved with a numerical fit to obtain the individual masses [|], [|, |j. It was recently shown 
that the mass precision obtained from this technique can be improved by subsequently per- 
forming combined fits to individual events, imposing both experiment end-point constraints 
and event E™ lss constraints ||]. 

When the number of kinematic end-point constraints provided by a given decay chain 
is insufficient to fully constrain individual sparticle masses alternative techniques must be 
employed. One possible approach involves solving simultaneously the mass-shell conditions 
obtained from several events containing the same decay chain This mass-relation 
method exploits the small widths of SUSY states, allowing the mass of each state appearing 
in the considered events to be assumed to be constant. 
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A second approach to this problem is to select events in which the same decay chain 
appears in both 'legs' of each selected event. In this case additional constraints are pro- 
vided by the components of the event Ex"" ss vector 1 , and again use can be made of the 
sparticle narrow-width approximation to equate the masses in the two legs. This permits 
the construction of further distributions with kinematic end-points related to the masses 
of sparticles present in the event. 

One example of a technique of this kind is the stransverse mass method || [jj. Con- 
sider two identical heavy SUSY states S± and 82 , decaying respectively to visible products 
v\ and V2 and identical lighter states a.\ and 02- If Pt(o;i), the transverse momentum 
vector of a\, were known then it would be possible to calculate t«t(<5i), the transverse 
mass of 5i, which is bounded from above by m(5). If Pt(cki) is known however then so is 
Pt(o ; 2) through application of the event Ex™ 55 constraints. Therefore my (£2) could also 
be calculated, a quantity which must also be less than m(5). Consequently the maximum 
value of tot(<5i) and mxi^) provides a variable with an end-point whose position measures 
m(5). Of course in reality we are not able to measure pt(cki) or Pt(«2) however the great 
insight of Ref. || was the realisation that if we can find a test value Pt(oi) which minimises 
this maximum transverse mass, we can be sure that the minimised-maximised transverse 
mass is also bounded from above by m{5). This 'minimax' transverse mass quantity is 
referred to as the 'stransverse mass' or Mt2- 

The development of the stransverse mass technique was particularly important because 
for the first time it allowed the measurement of masses of sparticles decaying through very 
short cascades, for instance cjr — ► qxi or I — > Zx?- Furthermore an analytical expression 
for Mt2 has recently been derived, valid in cases where the centre-of-mass (CoM) frame 
is at rest in the laboratory transverse plane ||, thus simplifying its use considerably. The 
technique inherits one draw-back from its use of the transverse masses of 5 decay products 
however, namely that it requires the use of m(a) as an input. Mt2 may therefore be 
described more correctly as an ensemble of variables, one for each assumed value for the 
unknown quantity m(a). The dependence of Mt2 on m(a) has been determined to be 
approximately m(5) — m(a) in specific cases [Q], however it would in general be preferable 
if the definition of the variable were independent of the unknown quantities to be measured. 
In that case the mass constraints obtained from an end-point fit would be uncorrelated with 
other measurements and hence could be used as input to a global mass fit. 

In this paper we will propose a very simple technique which seeks to address the same 
problem as the stransverse mass technique, but which approaches the problem from a 
different perspective. The new technique will allow a simple analytical combination of par- 
ticle/sparticle masses to be constrained in a precise and model- independent manner. Fur- 
thermore the technique will offer at least in principle the prospect of measuring individual 
particle masses, as postulated for the stransverse mass technique in Refs. ||, [11], O, |i~2| . 
The new technique will be applicable to any class of events in which heavy particles of 
unknown mass are pair-produced and decay to semi-invisible final states. 

1 We denote three- vector and two-vector quantities with bold case, while the corresponding magnitudes 
are denoted with standard case. Four-vector quantities are written in standard case. 
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The structure of the paper is as follows. Section |2| will describe the principles underlying 
the technique and investigate the properties of the new variable upon which it is based. 
Section || will illustrate application of the technique to the problem of constraining sparticle 
mass combinations with q~R — > qx± pair events at the LHC. Section |I| will outline extension 
of the technique to measurement of individual particle masses. Section || will conclude and 
discuss avenues for future work. 

2. Description of technique 
2.1 Background 

Consider 'symmetric' events in which identical cascade decay chains of the form 

5->av (2.1) 

occur in each leg i of the event. We shall refer to the initial particles produced in the hard 
scattering as Si. We shall further consider n step decay chains in each leg consisting of 
n — 1 decays, such that the (n — l) th decays produce invisible particles a,. The visible 
products of decays 1 to (n — 1) in each leg will be considered as single systems Vi of mass 
m{vi) and four-momentum p(vi). We shall assume that no invisible particles other than a 
are produced in the decay chains. The particles Si and on have common masses which are 
respectively m{S) and m(a). 

This parameterisation of the decay chains is quite general. The case n = 2 corresponds 
to SUSY chains such as qR — > qxi or / — > Ixi, with a identified as the LSP Xi- I n these 
cases m 2 (vi) « p 2 (vi) at the LHC. For longer SUSY chains we can choose the number 
of decays provided we can unambiguously identify the visible products of those decays. 
If n is equal to the total number of sparticles in the chain then a is again the LSP. For 
chains with n > 2 steps the distributions of invariant masses m(vi) can display kinematic 
end-points sensitive to analytical combinations of sparticle masses appearing in the chain 
Q. This information is used to constrain the individual masses in the end-point method 
but will be incidental to the technique described here. 



Consider now the use of iV symmetric events of the general form of Eqn. to mea- 
sure m{S) and m(a). This problem reduces to one of solving 6iV non-linear simultaneous 
equations, with each event providing the mass-shell conditions: 

[p(vi) +p(a 1 )} 2 = \p(v 2 ) +p(a 2 )} 2 = m 2 (S), 

\p( ai )] 2 = \p(a 2 )) 2 = m 2 (a\ (2.2) 

together with two E, T miss constraints: 

p x (ai)+ Px {a 2 ) = E™ ss , 

p y {a l )+p y {a 2 ) = E™ ss . (2.3) 

Each event contributes 2 unknown masses, which are common to all events, and 8 unknown 
on four-momentum components, which differ between events. The total number of unknown 
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parameters is therefore 8iV + 2 while the number of constraints is 6iV and so the system 
of equations is highly under-constrained. 

It may seem surprising at first that the above system of equations can be solved at all, 
however it should be noticed that we are not concerned with measuring the four-momenta 
p{a\) and piu-i) for all events, but rather with measuring only the common masses m(5) 
and m(a) using at least one event. Consequently we may set out to discard events in 
which the unknown masses depend on unknown four-momentum components. The problem 
therefore reduces to one of finding variables dependent only on the measurable quantities 
p(v\), p(v2) and Ex™ ss which identify events where the masses also depend only on those 
measurable quantities. This general approach is effectively that taken by kinematic end- 
point techniques, in which the variables identifying the events, such as mill) or Mt2, are 
also those which provide the mass measurment. This is also the approach which shall be 
taken here. 

2.2 Transverse momentum end-points 

One possible starting point for this problem was outlined in Ref. |^]. In an effective two- 
body decay process of the type considered above the magnitude of the three-momentum 
of the visible decay products in the rest frame of Si is given by 



\p(vi 



1 \J[m 2 (5) — m 2 (a) + m 2 (v j)] 2 — [2m(S)m(vi)] 2 
= 2 m{5) 

= \M h (2.4) 



which defines the 2-body mass parameter Mi. It will also be useful for the discussion 
which follows to define the equivalent quantity Mq for the special case where m(vi) = 0: 

-Mo - m2{5) -f (a) - (2.5) 
m(o) 

If Si has a small boost in the laboratory transverse frame, then the laboratory transverse 
momentum of Uj is of order Mi/2. This dependence of the momenta of visible decay 
products on the masses of heavy particles further up the decay chain is the reason that 
variables such as the 'effective mass' Q used in SUSY studies are sensitive to such masses. 

In principle we can improve on the use of ad hoc variables such as the effective mass 
however. In the rest frame of Si the magnitude of the momentum of Vi transverse to the 
beam direction 2 , pr(vi), is bounded from above by Mi/2 because 

PriVi) = ^-sin^i, (2.6) 

where ipi is the polar decay angle relative to the beam direction. Consequently if we could 
measure PT{vi) we could constrain the masses. 

Unfortunately however we are not able to measure priyi) directly - instead we measure 
the equivalent quantity in the laboratory frame: Px(vi). To proceed further we assume 



2 We denote quantities measured in the S1S2 CoM frame with primed variables and those measured in 
the rest frames of Si or <$2 with unprimed variables. 
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that the <5i<52 CoM frame is at rest in the laboratory transverse plane. This condition can 
be enforced by selecting events in which the net transverse momentum of the final state 
excluding the <5i and 82 decay products is small. In this case p' T {vi) is related to px{vi) by 
a proper Lorentz transformation in the transverse plane through the well-known relation: 

PtM = -^Tp2 \PT(vi) cos & ± PE( Vi )] 2 +p 2 T (vi) sin 2 (2.7) 

where (3 is the transverse boost factor (0 < (3 < 1), E(vi) is the energy of vi and fa is the 
angle in the rest frame of 8i between the boost direction and Pt(^)- For given V{ we know 
neither (5 nor fa and hence we are not able to reconstruct pr{vi)- Nevertheless we do know 
from conservation of momentum that in the 8182 CoM frame, and hence the laboratory 
transverse plane, the boost applied to V2 is equal and opposite to that applied to v\. 

To proceed further we shall attempt to find a quantity which can be calculated from 
the components of Pt(^i) and pt(^2) which remains unchanged if calculated with the 
corresponding components of p^C^i) and y>'t{ v 2)- If we could find such a quantity then we 
could use it to relate momenta measured in the 8182 CoM frame to those measured in the 
8\ and 82 rest frames and hence constrain Mi. 

2.3 Cotransverse mass and contransverse mass 

Consider first a system containing two particles v\ and i> 2 with masses m{v\) and m(v 2 ) 
measured in some frame F( ) to have four-momenta p{v\) and p(v2)- If both these particles 
are now measured in a different frame F^) it is well known that the mass obtained from 
p{v\) + p{v2) remains unchanged, i.e. the quantity 

m 2 ( Vl , v 2 ) = [E( Vl ) + E(v 2 )] 2 - [p(vi) + p(v 2 )] 2 

= m 2 {vi) + m 2 (v 2 ) + 2[E(v 1 )E(v 2 ) - p(vi) ■ p{v 2 )] (2.8) 

is invariant. Another way to interpret this is that when particles v\ and i>2 are subjected 
to co-linear boosts of equal magnitude m 2 (vi, v 2 ) is invariant. 

Now let us examine what happens when we start from one frame F( ), but boost 
particles v\ and V2 to different frames F^) and F( 2 ) respectively. These new frames are 
distinguished by the fact that their boosts are of equal magnitude but opposite direction in 
frame F/ ). In other words particles v\ and v 2 are subjected to contra-linear boosts of equal 
magnitude. Clearly m 2 (v\,V2) is no longer an invariant - this can be seen for instance by 
considering p(v\) = — p(«2) i n which case E(v\) + E{v2) increases with increasing while 
p(v±) + p(t>2) remains zero. 

Consider now a new quantity Mq equivalent to the invariant mass obtained from 
p{v\) + V{p(v2)) where V is the standard parity transformation operator: 

M&VUV2) = [E{ Vl ) + E{v 2 )) 2 - \p(vi)-p(v2)] 2 

= m 2 { Vl ) + m 2 (v 2 ) + 2[E(v 1 )E(v 2 ) + p(«i) ■ p(v 2 )}- (2.9) 

This quantity is invariant under the contra-linear boosts considered above. Denoting quan- 
tities measured in F( ) with primed variables, and those measured in F^) and F( 2 ) with 
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unprimed variables, and defining the x direction to be the boost direction, this can easily 
be demonstrated: 



M'S( Vl ,v 2 ) = [E'( Vl )+E\v 2 )) 2 - [p'M -p> 2 )] 2 

= 7 2 [E(vx) + Px {v x ) + E(v 2 ) - f3p x {v 2 )] 2 
-I 2 \Px{vi) + (3E{v x ) - Px {v 2 ) + (3E{v 2 )f 

~ \Py(vi) -Py(v 2 )] 2 ~ [Pz(vi) ~Pz(v 2 )] 2 

= 7 2 ([E( Vl ) + E{v 2 )] 2 + (3 2 [p x ( Vl ) - p x {v 2 )) 2 + 2p[E(vi) + £(«2)][ps(ui) - p x (v 2 ))) 
- 7 2 ((3 2 [E( Vl ) + E(v 2 )} 2 + \p x (vi) -p x (v 2 )) 2 + 2(3[E( Vl ) + E(v 2 )}[p x ( Vl ) - p x (v 2 ))) 

~ \Py(vi) ~Py(v 2 )] 2 ~ [Pz(vi) -Pz{v 2 )f 

= 7 2 {[E( V1 ) + E(v 2 )] 2 [l - f3 2 ] - \p x {v x ) - Px (v 2 )] 2 [l - P 2 ]) 

~ \Py(vi) -Py{v 2 )} 2 ~ [Pz(vi) ~Pz(v 2 )] 2 

= [Eiv^+Eiv^-lpiv^-p^)] 2 

= M%(v 1 ,v 2 ). (2.10) 

Since Mc(vi,v 2 ) is invariant under contra-linear boosts of equal magnitude its value can 
be calculated from the momenta and energies of v% and v 2 in any pair of frames and 
F( 2 ) related to F( ) by such boosts. For instance in the case considered above Fm\ could 
be identified with the 5\5 2 CoM frame and and F( 2 ) identified with the rest frames of 
5i and 5 2 , in which |p(i>i)| = M.i/2 and |p(^2)| = M. 2 /2. 



From a practical perspective the quantity Mc{v\,v 2 ) defined by Eqn. 2.9 is relevant 
only to cases where the 5i5 2 CoM frame is at rest in the laboratory frame, for instance in 
collisions at a lepton collider such as LEP or the ILC. At a hadron collider the scenario 
is more complicated. As discussed above, co-linear boosts in the laboratory transverse 
plane can be limited by selecting events in which the net transverse momentum of the 
final state excluding the 5\ and 5 2 decay products is small. There remains however a 
potentially large co-linear boost in the beam (z) direction caused by the differing proton 
momentum fractions of the colliding partons in the event initial state. Mc(vi,v 2 ) is not 
invariant under co-linear boosts of v\ and v 2 because V does not commute with proper 
Lorentz transformations. Consequently we must focus purely on quantities constructed 
from momentum components measured in the laboratory plane transverse to the beam 
direction. 

If v\ and v 2 were subjected to co-linear rather than contra-linear equal magnitude 
boosts in the laboratory transverse plane then a suitable invariant quantity to consider 
would be the transverse mass mxivi, v 2 ) [^4|, hereafter refered to as the cotransverse mass. 
mT(v±,v 2 ) is defined by: 

rriT(vi,v 2 ) = [E T (vi) + E T (v 2 )] 2 - [pt^i) + Vt{v 2 )] 2 

= m 2 ( Vl ) + m 2 (v 2 ) + 2[E T (y x )E T {v 2 ) - p T («i) • PtM], (2.11) 



-6- 



where 

E T (vi) = y/l%(vi)+m 2 (vi). (2.12) 

This quantity is useful because it is bounded from above by m(vi,v 2 ). When m{v\) = 
fn{v 2 ) = the following simplification can be made: 

m^(v 1 ,v 2 ) = 2p T (vi)p T (v 2 )(l - cos ^12), (2.13) 

where <pi 2 is the angle between v\ and v 2 in the transverse plane. This illustrates that 
events saturating the bound on the (co)transverse mass typically require that v\ and v 2 be 
back-to-back. 

In the case of contra-linear equal magnitude boosts considered above the equivalent 
quantity to the (co)transverse mass can be derived from Eqn. 2.9 : 

Mq T (vi,v 2 ) = [E T {vi) + E T (v 2 )] 2 - [p T (vi) -PtC^)] 2 

= m 2 («i) + m 2 (v 2 ) + 2[E T ( Vl )E T (v 2 ) + p T («i) • Pt{v 2 )]. (2.14) 

We shall refer to this quantity as the contransverse mass. This has the property that when 
m{v\) = m(v 2 ) = it reduces to 

Mq T (vi,v 2 ) = 2p T (vi)p T (v 2 )(l + cos0i 2 ), (2.15) 

where if Pt{v\) and pt(v 2 ) are measured in the laboratory transverse plane then <pi 2 is the 
angle between v% and v 2 in that plane. 

It is interesting to note at this point that when v\ and v 2 are massless and the <5i<52 CoM 
frame is at rest in the laboratory transverse plane the 'E T miSS vector can be represented 
under a change of basis involving Mct(vi,v 2 ) : 

E T miss = { _ Px ( Vl) -p^^-pyivx) -p y (v 2 )} - {pt(vi)-Mv2),Mct{v x ,V2)}. (2.16) 

In the new basis the first component can be interpreted as the contribution to E™ lss from 
Pt asymmetry, while the second, containing the geometric mean of pt{v\) and Pt(v 2 ), can 
be interpreted as the contribution from event topology. In this case E™ lss is given by: 



ET SS = \J\Pt{vi) -Pt(v 2 )} 2 + M* T ( Vl ,v 2 ). (2.17) 

The physical interpretation of the contransverse mass is more difficult than in the 
(co)transverse case. Mct(vi,v 2 ) does not represent the mass of a particle decaying to 
produce v\ and v 2 . Nevertheless we expect its distribution to display an end-point because 
it can in principle be calculated from the momenta of visible decay products measured 
in the rest frames of 5\ and 5 2 , and we know from Section |2.2| that these momenta are 
bounded from above by Mi/2. For instance if m(v±) = m(v 2 ) = then Mct{v\,v 2 ) takes 
a maximum value of Mo, i.e. 

= m2{5) ~f {a) . (2.18) 
m(b) 
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Interestingly this bound is saturated when v\ and V2 are co-linear, in contrast to the case 
for the (co)transverse mass. 

To summarise, we have now found a quantity bounded from above by an analytical 
combination of particle masses, and which can be calculated using momenta of visible decay 
products measured in the laboratory transverse plane. We shall now consider as a use-case 
the practical application of this variable to LHC data in order to measure SUSY particle 
masses. 



3. Example: q^qji events at the LHC 

To illustrate the application of the contransverse mass end-point technique to LHC data, 
a Monte Carlo simulation study was carried out aimed at measuring M™£ x for q~R pair 
production events where each q~R decays to a quark and a Xv Squark mass measurement 
in this channel using the stransverse mass method was first studied in Ref. |jl5| . The 
experimental signature of this process is the presence of events with exactly two jets and 



large E™ tss . In the context of the decay chain discussed in Section 2.1 the q~R plays the 
role of 8 and Xi that of a. We assume that the quark jet decay products are massless 



and hence Eqn. [2.18| allows us to measure an analytical combination of sparticle masses by 
measuring M^ x . 

A sample of 480k SUSY signal events equivalent to 10 fb _1 of data was generated 



from the SPSla benchmark mSUGRA model |L§ with HERWIG 6 . 5 |T§ |7| and passed to 
a generic LHC detector simulation |l^] modified to impose an 80% efficiency for electron 
identification, with mis-identified electrons being added to the list of jets if pr(e) > 10 
GeV. The ISASUGRA 7.69 RGE code |9| was used to calculate the input SUSY mass 
spectrum, giving m{qn) ~ 548 GeV and m(x?) = 96 GeV and hence M™f x = 531 GeV. A 
fully inclusive sample of SUSY events was generated in order to model SUSY backgrounds. 
Events were selected with the following requirements (with ji used to denote jet i): 

• rij e t = 2 for AR = 0.4 cone jets with pt{]) > 10 GeV and \r]\ < 5.0, 

• n iep = for isolated leptons (electrons or muons) with pt > 5 GeV (electrons) or 
Pt > 6 GeV (electrons), \r/\ < 2.5, minimum AR with nearest jet of 0.4 and maximum 
energy deposition of 10 GeV in a AR = 0.2 isolation cone, 

• min[p T (ji),p T (j 2 )] > 100 GeV, 

• EJf liss > 200 GeV, 

• in order to limit boosts of the q~Rq~R CoM frame in the laboratory transverse plane, 
measured pt of the j\j2 + EJp iss CoM frame in the laboratory transverse plane must 
satisfy 



\J[Px{h)+Px{h) + E™ iss } 2 + \Py{ji) + Py{j2) + E™^] 2 < 20GeV, (3.1) 
M CT > 200 GeV. 
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The hard Pr{j) an d E™ lss cuts additionally ensure that events easily pass typical LHC 
high level jet + EJp iss trigger criteria such as p T (j) > 70 GeV and EJp iss > 70 GeV @. 

After application of these cuts many Standard Model (SM) backgrounds are heavily 
suppressed: 

• QCD jet backgrounds, while possessing a very large cross-section, are suppressed by 
the Mct cut which rejects events with back-to-back jets. Jet energy mis-measurement 



mainly generates E™ tss through the first term in Eqn. 2.17 and the effect on Mct 
is smaller. In order to pass the Mqt cut at least one high px jet must be com- 
pletely missed by the detector. The Mct cut is strongly correlated with the D n7T 
variable used at the Tevatron to separate SUSY signal from QCD backgrounds in 



multijet+£ , ™ s,s searches [21]. The fast detector simulation used in this study is not 
expected to model the catastrophic loss of jets accurately, but we do not expect this 
background to be dominant even when using a more realistic simulation. In particu- 
lar such events can in principle be removed with 'event cleaning' cuts, for instance by 
reconstructing jets from charged particle tracks. Consequently QCD jet backgrounds 
are not considered further here. 

• Hadronic or semi-leptonic ti backgrounds are suppressed by the jet multiplicity cuts, 
while fully leptonic events in which both leptons are lost inside the jets (the worst 
case scenario kinematically) possess Mct values less than m(t) ~ 172 GeV, which is 
the value expected for top quarks decaying to a neutrino plus massless visible decay 
products. Such events therefore fail the Mct cut. 

• W + 1 jet backgrounds in which the W decays to a hadronic r or electron faking a jet 
can mimic 2-jet events. Events with large Mct typically possess two co-linear jets. 
When the lepton is emitted co-linearly with the initial jet its transverse momentum 
is given by 

PT = ±m(Wh(l-P), (3.2) 

where (3 is the boost of the W in the transverse plane. The maximum value of Mct 
generated by this configuration is obtained in the limit /3 — » 1, when M^ x = m(W). 
Such events therefore also fail the Mct cut and are not considered further here. 

The remaining SM backgrounds are dominated by Z(— > vv) + 2 jets and W{— > Iv) + 
2 jets events, where in the latter case the lepton momentum is anti-parallel to the W 
momentum and is 'red-shifted' such that its magnitude is below the lepton identification 
Pt threshold. These backgrounds were modelled with ALPGEN |z| coupled to HERWIG 6.5. 
In order to add realism to the analysis the backgrounds were estimated using data-driven 
techniques applied to the Monte Carlo 'data'. The Z(— > vv) + 2 jets background was 
estimated by selecting Z{— > It) + 2 jets events with similar cuts to those listed above, but 
replacing the lepton veto and EJp tss requirements with a requirement for two opposite-sign 
same-flavour leptons with \m{lt) - m(Z)\ < 10 GeV and |p T (Zi) + B T miss \ > 200 GeV. 
The W(— > Iv) + 2 jets background was estimated by selecting W(— > Iv) + 2 jets events in 
which the lepton was boosted and the neutrino de-boosted such that pr(t) > 200 GeV and 
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Figure 1: Distributions of Mqt values. The left-hand figure shows the cumulative 'data' dis- 
tributions summing SPSla SUSY events (light /yellow histogram), Z{— * vv) + 2 jets background 
events (medium/green), W(-^ Iv) + 2 jets background events (dark/blue) and diboson, top-pair and 
single-top background events (magenta/medium-dark). Data-points indicate the result of the data- 
driven estimate described in the text. The data-points in the figure on the right show the result of 
subtracting the data-driven estimate from the 'data' distribution, with the light/yellow histogram 
representing the SUSY distribution with no SM background added. The dark/red histogram shows 
the contribution from non-g^g^ SUSY events. The result of a simple linear end-point fit to the 
data-points is shown. 



j^jntss ^ -j^q Q e y E ac h data-driven estimate was normalised separately to the respective 
Monte Carlo background pt{]) distribution below the expected SUSY signal region. In 
practice the relative normalisation of the W(— > Iv) + 2 jets estimate could be obtained 
from data, for instance with a fit to the lepton pt spectrum in W(-^ Iv) + 2 jets events. 

SUSY backgrounds to qR(jR events arise primarily from processes in which at least one 
qi decays through a chain producing multiple invisible final state particles. One possible 
example involves sneutrinos decaying to neutrinos and Xi- I n these cases the mass of each 
SUSY state produced in association with the jet in the decay of each qi is greater than 
that of the Xi produced in the decay of <jr and consequently these events possess M™£ x 
values below those for q^qR. At parton level if SUSY background events are to exceed 
the expected end-point, assuming correct assignment of decay products to SUSY decay 
chains, then the mass of the initially produced sparticles must be greater than m(qR). The 
main candidate for this is g pair production in which each gluino decays to co-linear jets in 
association with a x.i- This process should generate a Mqt distribution with an endpoint 
at Mqj? = 597 GeV for a g mass of 612 GeV at SPSla. 

The Mct distribution for events satisfying the selection cuts is shown in Fig. |](left) 
indicating an excess of events at large Mqt values due to SUSY processes. As expected the 
contribution from ti events (modeled with HERWIG 6.5) is small, as are the contributions 
from WW, WZ, ZZ and single-top production (also modeled with HERWIG 6.5). The 



- 10 - 




200 



300 



400 



500 



600 



700 800 

M rT (GeV) 



Figure 2: Distributions of Mct values for SUSY signal events for ten different values of the cut 
on the measured pt of the jij2 + E™ lss CoM frame (Eqn. 3d), ranging from 200 GeV (top) to 
20 GeV (bottom) in 20 GeV steps. In contrast to Fig. [I] the jet multiplicity cut has been relaxed 
to require at least two jets and a logarithmic y-axis has been used to aid comparison of shapes of 
distributions. 



data-points in Fig. |T](right) represent the same distribution after subtracting the data- 
driven background estimate. As expected, a prominent end-point feature is visible at 
around 530 GeV. A simple linear fit to the endpoint determines its position to be 550 ± 53 
GeV (10% uncertainty). Use of a more sophisticated fitting function would undoubtedly 
improve this precision significantly. There is also some evidence in Fig. |l](right) for a small 
excess of events beyond the expected <jr(Jr end-point. Examination of the Monte Carlo 
truth record indicates that the large Mqt values of these events originate either from jet 
mis-measurement in qnqn events or from both jets originating from the same SUSY decay 
chain in non-quqn events. No gg events were observed to contribute to this region for this 
SPSla SUSY model. 

In order to study the dependence of the shape and position of the Mct end-point in 
Fig. m on the cut on the measured pt of the j'ij'2 + E™ 1SS CoM frame (Eqn. |3.1| ), the Mct 
distributions of SUSY signal events passing progressively harder pt cuts were generated. 
Due to the strong correlation between the di-jet multiplicity cut and the pt cut the former 
cut was relaxed to require at least two jets, with the two hardest jets being used to calculate 
Mct- The resulting distributions are plotted in Fig. [2| for ten different values of the pt cut 
ranging from 200 GeV (top) to 20 GeV (bottom). The effect of the cut on the SUSY signal 
is to sharpen the end-point at the expense of statistics. This sharpening of the end-point is 
caused both by limitation of the transverse boost of the ququ CoM frame and by rejection 
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of high multiplicity non-q/ji^R events in which two jets from the same SUSY decay chain 
are selected to calculate Mqt- It should also be noted that a harder px cut rejects more 
SM background events, especially QCD multijet events. 



4. Extension to measurement of individual particle masses 

So far we have shown that we can obtain an end-point in the distribution of a quantity 
calculated from visible decay product transverse momenta which depends on an analytical 
combination of masses. In principle it is possible to use the position of this end-point to 
measure the individual masses m{5) and m(a). The key to this is recognising that A4i 
depends upon both the unknown masses m{5) and m{a) and the visible masses m(vi). If 
one requires that m{v\) = m(v2) = m(v) then it can be shown from Eqns. 2A and [2.14 
that for given m(v) the bound on Mct is given by: 

M ct =-^m\v) + M Q . (4.1) 

Consequently if events could be found in which m{v\) and m{v2) were non-zero and 
equal, for instance by accurately combining products from several decays in a multi-step 
chain, then the position of M^ x would depend linearly on m 2 (v) with a gradient of 
l/m(5) and intercepting the ordinate at .Mo- The gradient of a linear fit to this bound 
therefore measures m{5) independently of m(a), while the intercept then allows m(a) to 
be constrained. As an aside it is interesting to note that in this case Mqt is also bounded 
from below by: 

M^ n = 2 V / ^2(V). (4.2) 

and for exclusive decay chains m 2 (v) is bounded from above by a separate analytical com- 
bination of masses identical to that used by the end-point method discussed in Section [l]. 

The above technique should work in principle however it is likely to be very difficult 
to implement in practice. The first difficulty is connected with unambiguously associating 
decay products with SUSY decay chains. One possible approach would involve focusing 
on specific exclusive decay chains, using the values of invariant masses of combinations of 
decay products to associate decay products to chains [Q]. Unfortunately however the low 
acceptance of such exclusive selections is likely to prevent successful application of this 
technique before significant quantities of data have been acquired. 

A second approach involves inclusive selection of SUSY events with multiple visible 
decay products, and use of a kinematic algorithm to approximately associate decay products 
to chains. Fig. ||(left) shows the result of an attempt at applying this approach to SPSla 
events with ISR turned off. Here decay products have been associated to chains by requiring 
that max[m 2 (wi), m 2 (f2)] is minimised. Events have been selected by requiring four jets 
and no leptons, and the mass-squared equality requirement mentioned above has been 
imposed by requiring that the asymmetry in m?(vi) and m 2 (t>2) is less than 10%. An 
additional cut requiring the rapidity difference between v\ and t>2 to be greater than 1.0 
has also been applied to reduce combinatorics. 
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Figure 3: Distribution in the Mct — rn 2 {v) plane of SPSla events with ISR turned off (left) and 
on (right). The top (straight) line represents the expected dependence from Eqn. 4.1 of M™^ x on 
m 2 (v) for events containing g pair production. The middle (straight) line represents the equivalent 
dependence for events containing q pair production. The bottom line represents the expected 
dependence of M™™ on rn 2 (v) given by Eqn. |4.2|. 



Fig. H|(left) shows that with these cuts events generally lie below the expected upper 
bounds on Mct for g decays (top line) or q decays (middle line), although some combi- 
natorial contamination is visible above the g bound. The lower bound on Mct given by 
Eqn. is prominent (bottom line). Fig. ||(right), obtained with SPSla events with ISR 
turned on, illustrates the further difficulty of using this approach however. The inclusion of 
ISR jets in v\ and/or V2 can artificially lower m 2 (v) and hence generate false configurations 
which strongly violate the upper bounds on Mct- A related effect was noted previously 
in connection with the stransverse mass technique in Ref. pi. Clearly much more work is 
needed before this technique can be used practically to measure accurately independent 
particle masses. 



5. Conclusions and directions for future work 

This paper has shown that by constructing a kinematic quantity invariant under contra- 
linear equal magnitude boosts in the laboratory transverse plane a simple analytical com- 
bination of the masses of pair-produced particles and their invisible decay products can 
be constrained at hadron colliders such as the LHC. It was shown that in principle these 
techniques may be used to measure the masses of such particles independently, although 
in practice this seems to be very difficult. 

The study described in this paper suggests several directions for future work. These 
include: 

• The experimental simulation study of Section || should be repeated with more realistic 



-13- 



full experiment-specific simulation of all Standard Model and SUSY backgrounds to 
demonstrate conclusively the feasibility of these techniques when applied to q~Rq~R 
events. The feasibility of I mass measurement with 11 events should also be studied. 

• Further work assessing the feasibility of measuring independent particle masses using 
the technique outlined in Section [|] is required, focusing in particular on optimising 
the experimental assignment of decay products to decay chains, and rejection of ISR 
jets. 
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